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The transient electric current of surfactants dissolved in a nonpolar solvent is investigated both experimentally 
and theoretically in the parallel-plate geometry. Due to a low concentration of free charges the cell can 
be completely polarized by an external voltage of several volts. In this state, all the charged micelles are 
compacted against the electrodes. After the voltage is set to zero the reverse current features a sharp discharge 
spike and a broad peak. This shape and its variation with the compacting voltage are reproduced in a one- 
dimensional drift-diffusion model. The model reveals the broad peak is formed by a competition between 
an increasing number of charges drifting back to the middle of the cell and a decreasing electric field that 
drives the motion. After complete polarization is achieved, the shape of the peak stops evolving with further 
increase of the compacting voltage. The spike-peak separation time grows logarithmically with the charge 
content in the bulk. The time peak is a useful measure of the micelle mobility. Time integration of the peak 
yields the total charge in the system. By measuring its variation with temperature, the activation energy of 
bulk charge generation has been found to be 0.126 eV. 
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I. INTRODUCTION 

Electrical conduction in nonpolar fluids is of interest 
in relation to electrophoretic displays that use nonpolar- 
based solutions as the working medium.— ~— When an ap- 
propriate surfactant is added to a nonpolar solvent its 
molecules form inverse micelles: globular aggregates with 
hydrophobic surfaces and hydrophilic interiors X The mi- 
celles are believed to be responsible for the solution's 
nonzero electrical conductivity by stabilizing positive and 
negative charges and keeping them from recombining, 5 
In addition, the micelles have been suggested to charge 
colloidal particles by ionizing their surface groups and 
preferentially adsorbing on the surfaces. 6 

A useful way to gain insight into the microscopic prop- 
erties of inverse micelles is through a detailed analysis of 
the transient current curves^— Typically, in such an ex- 
periment a single-surfactant solution is placed between 
two parallel electrodes, although more complex geome- 
tries have been usedJ^ The voltage is stepped from zero 
to a finite value and then either back to zero or to some 
other value. The current recorded in the external cir- 
cuit is an integral measure of the complex movement of 
charges inside the cell. By assuming the solution to be a 
symmetric electrolyte, the concentration of charged mi- 
celles, their mobility, diffusion coefficient, and the Stokes 
diameter can be inferred. 

Most of the previous work has focused on the some- 
times very detailed quantitative analysis of the forward 
transient current] 11 ' 16 However, in the early studies it 
was noticed that the reverse transients often exhibit a 
nonmonotonic behavior. 2,9 Although the reverse peak 
has been observed in subsequent studies) 13 ' 15 ' 18 it has 



not been analyzed or explained in sufficient detail. The 
present paper aims to contribute to this area. 



II. EXPERIMENT 

Poly-isobuthylene succinimide (OLOA 11000, Shevron 
Oronitc Company LLC, 72% active component) was dis- 
solved in the hydrocarbon fluid Isopar® M (Exxon Mo- 
bil Chemical, dielectric constant e — 2.0) by sonication at 
different concentrations. No attempt was made to dry the 
components and the water content was not controlled. 
Parallel-plate cells were made from pieces of polyethy- 
lene terephthalate (PET) coated with conducting indium 
tin oxide (Sheldahl, Inc.). One side of the cell was photo- 
patterned with a square grid of SU8 spacer walls whereas 
the other side remained flat. The width of the walls was 
10 /jm and the grid period was 230 fim. The height of 
the walls (10 or 30 /zm) defined the cell thickness. The 
active device area was about 10 cm 2 . The voltage pro- 
file was provided by a DC power supply (Agilent Tech- 
nologies). The current was recorded by a low- noise cur- 
rent preamplifier (Stanford Research Systems) at a sam- 
pling rate of 1000 Hz. Both instruments were controlled 
by a customized Lab View application (National Instru- 
ments Corp.). To suppress the Faradaic processes, the 
measurements were done at below-room temperatures of 
5°C, 10°C, and 15°C. The temperature was controlled by 
a convective environmental chamber (Autronic-Melchers 
GmbH) with accuracy of 0.1°C. 



III. OBSERVATIONS 



Transient currents of a 10 /im-thick 0.5 wt.% OLOA- 
11000 solution in Isopar M are shown in Fig. QJa). The 
Electronic mail: pavcl.kornilovich@hp.com voltage between two parallel electrodes was stepped to a 

Electronic mail: |yoocharn jeon@hp.com| nonzero value between 0.25 and 10.0 V for 15 s, and then 
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FIG. 1. Transient currents of poly-isobuthylene succinimide in Isopar M in a parallel-plate cell. The voltage is stepped from 
zero to V at seconds and then back to zero at 15 seconds, (a) The overall shape of the transient for 0.5 wt %. The cell gap 
is d = 10 /urn and the temperature is T = 10°C. The step voltages are 0.25, 0.5, 0.75, 1, 1.5, 2, 3, 4, 6, 8, and 10 V. (b) A 
magnified view of the reverse peak region. Note the saturation of the shape after 3 V. (c) The cell gap dependence for 0.5 wt 
%. V = 10 volts, (d) The V = 8.0 volts peaks for three temperatures of 5°C, 10°C, and 15°C. 



brought back to V for another 14 seconds. As observed 
by other investigator a 7 ' 9 ' 13 ' 15 i 18 the currents possess the 
following characteristic features: a fast decay in the for- 
ward direction, then a sharp spike followed by a broad 
peak in the reverse direction. Initially, the forward decay 
is fast but then it slows down. Phenomenologically, the 
forward transient curve is well described by an RC circuit 
with a decreasing number of charge carriers (increasing 
R). This issue is not considered in the present work. 
The forward current continues to evolve with increasing 
voltage: it grows larger in magnitude and shorter in du- 
ration. In contrast, the reverse current quickly saturates 
with voltage. The V = 2.0 and 10.0 V curves are nearly 
identical. Figure QJb) shows a magnified view of the re- 
verse peak area. The peak develops between 0.5 and 0.75 
V and saturates by V = 3 V. This suggests 3 V cause 
complete polarization of the cell: all available charges are 
compacted near the electrodes and the charge cloud no 
longer changes with a further increase in voltage. 



Variation of the transient current with the cell thick- 
ness has been investigated. The reverse current for two 
thicknesses 10 fim and 30 /mi are displayed in Fig. Hfc). 
As expected intuitively, in the thicker cell, the reverse 
peak develops more slowly. The ratio of the peak times, 
16.5/1.65 = 10.0, is close to the ratio of the cell thick- 
nesses squared, 3.0 2 = 9.0. This scaling is discussed in 
the following sections. Variation of the current with tem- 
perature has also been investigated. Figure LTId) com- 
pares the reverse transients for the three temperatures of 
5°C, 10°C, and 15°C. As the temperature decreases the 
peak shifts to longer times and reduces in height. 



IV. MODEL 

To develop an understanding of the observed be- 
havior, a one-dimensional drift-diffusion (Nernst-Planck- 
Poisson) model has been applied. The model is formu- 
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lated here in a dimensionless formjii Given the distance 
between two parallel electrodes d, the diffusion coeffi- 
cient of charged micelles D (both positive and negative), 
and the physical distance x and time t, the dimension- 
less distance and time are introduced as s = x/d and 
t = (D/d 2 )t. The electrostatic potential <f> is measured 
in units of the thermal voltage: -0 = {e/ksT) <j>. The 
micelles are assumed to be single-charged. The concen- 
tration of micelles n±(x,t) is measured in units of the 
equilibrium initial concentration of one sign, c±(s,t) — 
n±/no- The micelle mobility /i is related to the diffusion 
coefficient by the Einstein relation D = ksT/i. In such 
units, the micelle flux is given by a sum of the drift and 
diffusion contributions 



3± 



dc± _ dtfj 

OS OS 



(1) 



In the simplest version of the model, the dissociation- 
recombination processes in the bulk^ ' 19 i 20 and the Fara- 
daic processes at the electrodes^! are neglected. (Both 
are small for the experimental times and temperatures 
studied in this work.) As a result, the fluxes satisfy the 
local conservation laws 



dc± 



ds 



= 0. 



and the blocking boundary conditions 
i±(0,r) =j ± (l, r )=0 



(2) 



(3) 



These conditions also imply global conservation of posi- 
tive and negative charges. The electrostatic potential is 
governed by the Poisson equation 



^=-X(c + -c_), 



(4) 



and boundary conditions 

#),r)=0; ^(l,r) = ^ = u(r). (5) 

k B T 

The step-function boundary conditions arc defined as 



0, r < Tl 
u(t) — \ u , Ti < r <t 2 
0, r > r 2 



(6) 



In Eq. (HJ), K — e 2 d 2 no/ 'eeoksT is a dimensionless 
parameter which is an equivalent measure of the to- 
tal charge concentration. It can also be expressed as 
K = d 2 /2L 2 J , where Ljj is the Debye length of the undis- 
turbed solution. Another parameter is the magnitude of 
the voltage step Uo- The model is fully characterized by 
just two dimensionless parameters, K and uq. As such, it 
allows for complete numerical investigation. The initial 
conditions are uniform concentrations and a zero electro- 
static potential: 
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FIG. 2. Numerical solution of equations H])-([7|) for different 
values of K and Mo- The voltage is stepped to uo at n =0.1 
and back to zero at T2 = 0.6. In all of the panels, Mo increases 
from 5 to 40 in steps of 5. Notice the saturation of the reverse 
peak with uq. Legend numbers apply to all panels. 



Once ip(s,T) is known, the external current (in amperes) 
is computed as the time derivative of the electric field at 
the electrode (as follows from the Gauss law): 



I(t) 



ee ADk B T d 2 ip{s,r) 



:d 3 



ds di 



s=0 



c+(s,0) = c_(s,0) = 1 



VM)=0. (7) 



where A is the electrode area. 

The Cauchy problem (HJ-® has been solved numeri- 
cally using MATLAB (The MathWorks, Inc). The solu- 
tion for several values of K and uq is shown in Fig. [21 
The model transient possesses the same qualitative fea- 
tures as the experimental one: a fast monotonic decay in 
the forward direction, then a sharp spike followed by a 
broad peak in the reverse direction. The reverse current 
saturates with uq- The shape of the reverse peak slowly 
evolves with increasing K. 

For the following discussion it is instructive to rewrite 
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formula (JSJ) to make the relationship between the exter- 
nal current and the distribution of charges inside the cell 
more explicit. Integrating by parts and applying the 
Poisson equation (UJ) one obtains for the potential dif- 
ference between the electrodes 



u{t) = 



ds 



K 



(c_|_ — c_) ds . 

(9) 



Note that due to global electroneutrality the electric 
fields near the two electrodes are always equal, or 
[dip(s,T)/ds] s = = [dip(s,T)/ds] s= x. Substituting © 
into (El one obtains 



I(r) du{r) 



J dr K Jo 3 \ i)r iJt 
The dimensional form of the last formula reads 



f dc + (s,r) <9c_(s, r) 



ds . 

(10) 



x [n + (x,t) — n_ (x,t)]dx , 
( U ) 

where the dot means a partial derivative with respect 
to the real time t. The external current comprises two 
different contributions. The first term is charging and 
discharging of the geometric capacitor. (The coefficient 
by V is recognized as the geometric capacitance of the 
cell C g .) The second term is an additional current as- 
sociated with the movement of ions inside the cell. The 
first process is typically much faster than the second, and 
results in sharp charging/discharging initial spikes. 

Applying the conservation laws (|2]l. integrating by 
parts, and making use of the blocking boundary condi- 
tions (B]), formula (fTU|) can be rewritten as follows 

7(r) _ du(r) _ R /' 1 [j+(SjT) „ J _ (SjT)]ds . (12) 



dr 



Thus the long-term part of the external current is simply 
an integral of the ionic flux over the cell thickness. 



V. ORIGIN OF THE REVERSE PEAK 

In addition to the external current, the numerical solu- 
tion provides time-dependent distributions of the charges 
and the fields inside the cell, some of which are shown in 
Fig. [3l The origin of the reverse peak can be analyzed in 
sufficient detail, which is now described. 

During the forward phase the external electric field 
polarizes the cell by compacting the positive and nega- 
tive ions against their respective electrodes. The equilib- 
rium concentration profile is determined by the balance 
between the electrophoretic drift and diffusion when it 
reaches an equilibrium, cf. Fig. EJa), r = 0.1750. As 
the field is turned off, the electrodes are brought to the 
same potential on a time scale much shorter than the 
characteristic diffusion time of the ions. A new potential 
distribution is established almost instantaneously after 
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FIG. 3. Spatial and temporal dependence of physical quanti- 
ties for K = 800 and uo = 50. (a) Concentration of positive 
ions, (b) Electrostatic field, (c) Electrostatic potential. All 
quantities are dimensionless. The legend numbers are dimen- 
sionless times t. The legends are common for all three panels. 
The field is switched on at r = 0.1 and off at r = 0.6. The 
concentration of the negative ions is equal to that of the pos- 
itive ions reflected around the cell center s — 1/2. 



the switching at r = 0.6, cf. Fig.Efc), r = 0.6015. This 
fast geometric capacitance discharge produces the sharp 
initial spike. 

Just after the switching most ions are still being pushed 
against the electrodes by the local electric field. How- 
ever, the electrostatic potential of both electrodes is now 
zero. By symmetry, the potential in the cell center is 
also zero. (In a more general case of the asymmetric 
electrolyte or cell, the zero-potential point will not be 
in the center of the cell. This fact does not change the 
following arguments.) Therefore the potential assumes 
a local maximum somewhere near the negative electrode 
(and a local minimum near the positive electrode). This 
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zero-field point must be within the charge cloud since by 
Poisson equation only uncompensated charge can change 
the curvature of the potential. As a result, to the right of 
the maximum point there is some positive charge, which 
is of the order of the geometric capacitance times the 
minimal voltage of total polarization. (The estimate fol- 
lows from the Gauss law.) Since the electric field is also 
positive there, the charge-field product creates a flux of 
positive ions toward the center of the cell. A similar drift 
of the negative ions away from the positive electrode is 
taking place at the opposite side of the cell. Closer to 
the electrodes, the drift and diffusion fluxes continue to 
effectively compensate each other. 

Initially, the integral ionic flux is small and, according 
to Eq. (|T^|) . the external current is also small. However, 
the external current is a flow of the electrode charge, and 
as a result the latter decreases. The reduction of the 
electrode charge moves the zero-field point deeper into 
the charge cloud thereby releasing more ions for drift- 
diffusion toward the bulk. The electric field increases in- 
side the cell for the same reason. The integral ionic flux 
grows and so does the external current. Thus the current 
(i.e. time decrement of the electrode charge) is linearly 
related to the electrode charge itself. As a result, the cur- 
rent grows approximately exponentially, at least for short 
times when the positive and negative charge clouds do 
not overlap. As the clouds expand, positive and negative 
ions begin to mix in the middle of the cell reducing the 
space-charge and the electric field. The charge distribu- 
tion becomes more uniform which dampens the internal 
fluxes and the external current. The entire process can be 
viewed as a competition between the increasing number 
of participating ions and decreasing driving fields. In the 
early stages of the reversal the increase in the ion number 
dominates, and the overall current increases with time. 
Later, the ion number levels off but the driving electric 
field subsides as does the current. 

The variation of the reverse current with the charge 
parameter K oc d 2 no is now discussed. As can be seen 
in Fig. [21 larger Ks require larger voltages to reach cur- 
rent saturation. In addition, the peak maximum shifts 
to longer times. The voltage increase can be understood 
from the following considerations. The reverse current 
saturates when the cell is fully polarized which occurs 
when CdiV > Q. Here Cdi oc ^/Tiq exp (u p /A) is the 
Gouy-Chapman double-layer capacitance and Q oc nod. 
At low charge densities, the polarization voltage scales as 
Up oc dy/n~o oc \J~K. At moderate and high densities the 
exponential enhancement of the capacitance is most im- 
portant, and the scaling crosses over to u p oc ln(A"). It is 
essential that the polarization voltage scales sublinearly 
with the density. The above estimates do not take into 
account the significant effects of the limited charge con- 
tent in a finite cell on the double layer capacitance and 
polarization voltage. This issue was thoroughly analyzed 
by Verschueren et al.— 

As mentioned above, the reverse drift starts with a 
small amount of charge and then accelerates exponen- 
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FIG. 4. The saturated reverse current (uo = 50) at different 
values of the charge density parameter K. The reverse cur- 
rent is completely extinguished by r = 0.2. The curves have 
negative skewness whose magnitude increases with K. Notice 
a relatively weak dependence of the peak position on K. 



tially by pulling ions from the cloud near the electrodes. 
The process continues until the charge "stocks" are de- 
pleted. The depletion time approximately corresponds to 
the peak time in the reverse current. Therefore the peak 
time should scale logarithmically with the total amount 
of charge in the system, t p oc ln(A"). This scaling is con- 
firmed by numerical calculations, as shown in Fig.|4] No- 
tice the envelope of the curves' leading edges is roughly 
an exponential. 



VI. MICELLE PROPERTIES 

The broad peak is the most conspicuous feature of the 
reverse transient current in the high-voltage regime. In 
this section it will be shown how useful physical infor- 
mation can be extracted from the time and shape of the 
peak. 

Diffusivity. As can be seen in Fig. 21 the weak logarith- 
mic dependence places the peak time within the interval 
t, p = (0.10 ± 0.05)(gP/L>) for a large spread of the den- 
sity parameter K (a factor of 30) . This suggests a simple 
method of estimating the ion's diffusivity: measure the 
reverse peak time t p in seconds, then D ss 0.1(d 2 /t p ) 
with the accuracy of 50%. For example, for the data 
of Fig. md), the peak times are 2.00, 1.65, and 1.40 
sec for T = 5°C, 10°C, and 15°C. Using d = 10.0 //m, 
one obtains the respective diffusivities of OLOA11000 in 
Isopar M at these temperatures as D sa 5.0, 6.1, and 7.1 
/im 2 /sec. By the same argument, the peak time should 
scale approximately quadratically with the cell thickness 
d, cf. Fig.LUc). 

A more accurate method involves taking into account 
the changes of the peak shape with K. A convenient mea- 
sure of the current asymmetry is its skewness. Figure [S] 
shows the peak time at saturation versus skewness, as 
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FIG. 5. The dimensionless peak time of the saturated (i.e. at 
large uq) reverse current as a function of the skewness. The 
time is counted from the moment the voltage is set to zero. 



calculated from the model data of Fig. @] To avoid ambi- 
guity in determining the skewness the data set is limited 
to a symmetric interval around the peak time. Once the 
skewness is calculated for an experimental reverse current 
the dimensionless peak time t p can be found in Fig. [5] 
Then the micelle diffusivity is D = T p (d 2 /t p ). The skew- 
ness of the experimental data of Fig. [1] are given in the 
fourth column of Table HI The next column contains the 
respective dimensionless peak times t p determined from 
Fig. EJ They are 11-13 % larger than the "back-of-the- 
envelope" value 0.1, which corrects the diffusivities to 
their final values given in the sixth column of the table. 

Micelles Stokes mobility and Stokes diameter. Once 
the diffusivity is known the Stokes mobility fj,$ — 
D/kBT a and Stokes diameter 



1 



as = 



3irrj fig Snr] D ' 



(13) 



can be easily determined. Here T a is the absolute temper- 
ature in K, and r\ is the viscosity of the host fluid. The 
viscosity of Isopar M has been measured for — 25° C < 
T < +50°C. Note that the Stokes mobility differs from 
the more familiar electrophoretic mobility by the elemen- 
tary charge factor, [i ep — e us- The values of «s, a s an d 
rj for three temperatures are given in Table Q] Both us 
and rj show significant temperature dependence while the 
Stokes diameter does not. This suggests that the bulk of 
the temperature variation in the micelle mobility is due 
to the corresponding change of the viscosity of the host 
fluid. A similar conclusion was reached by Novotnyj^ The 
overall estimate of the micelle diameter is 12.2(5) nm. 

A similar analysis of the d — 30 «m sample at T = 
10°C resulted in a slightly larger mobility of 1.87 • 10~ 12 
s/kg and a slightly smaller micelle diameter of 11.8 nm. 
Note that the reverse transient current, cf. Fig. Hfc), 
is more skewed than the d = 10 /jm curves. The cor- 
responding skewness parameter, -0.76, is about twice as 
large. 



Equilibrium charge density. The saturated reverse 
transient curve allows an independent determination of 
the total amount of charge in the system. The Faradaic 
processes are suppressed at low temperatures and the 
bulk charge generation is too slow to affect the overall 
charge dynamics on the time scale of seconds^ Under 
these assumptions, the total amount of the charge of each 
polarity is conserved. This can be determined from the 
experimental transient curves. The argument is based on 
Eq. (jlip . By integrating the second term from the mo- 
ment just after the voltage turn off t\ (i.e., just after the 
discharge spike) to infinity, one obtains 



I(t) dt 



eA 
~d 



x [n + (x, <i) — n_{x, t\)} dx . (14) 



The t — oo contribution vanishes since n + (x, oo) = 
oo) = no, the micelles are fully mixed, and the 
integrand is identically zero. The volumetric densities 
n±(x, t\) are those just after the voltage turn off. In 
the ideal case of complete compaction, the densities are 
delta functions: n + (x,t\) = (nod)5(x), n-{x,t\) — 
(nod) 5(x—d). Performing the integration, the right hand 
side of Eq. ([TO becomes — e(Ad)nQ. Since Ad is the cell 
volume, the last expression is the total negative charge 
of the cell. In practice, the micelles are never compacted 
into a delta function; one can only speak about differ- 
ent degrees of approaching the limit. (In the experiment 
presented here, the compaction is believed to be almost 
complete, within a few percent of the limit.) The current 
integral provides the lower limit for the total charge: 



|Q+| = |Q-|> 



lit) dt 



(15) 



Once Q± are known, the equilibrium concentrations of 
charged micelles are found by dividing by the cell volume. 
The values of Q and no are given in the last two columns 
of Table |U The order of magnitude is no ~ 10~ 8 mol/1 
~ 10 -5 mol/m 3 , which is similar to the values reported 
in the literature^ 

Activation energy. The temperature dependence of the 
equilibrium charge density allows us to estimate the acti- 
vation energy E a of the charge generation process, which 
is of fundamental importance in dielectric fluids. Fitting 
the no(T) values from Table |H the activation energy is 
determined as E a — 0.126 eV = 2.91 kcal/mol. This value 
is close to the one reported by Novotnyii for Aerosol OT 
surfactant (0.140 eV), which suggest that the ionization 
mechanisms in the two systems might be similar. 



VII. 



SUMMARY 



In summary, transient currents in a weakly conduct- 
ing solution of poly-isobuthylene succinimide in Isopar 
M have been investigated. The limited charge content 
allows the complete polarization of the cell at low and 
medium voltages. Under these conditions, the bulk of 
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TABLE I. Properties of inverse micelles of OLOA11000 in Isopar M extracted from the reverse transient peaks, d is the cell 
thickness. T is the experiment temperature. t p is the time of the reverse current peak counted from the voltage turn off. The 
fourth column is the skewness calculated from the data interval symmetric around the peak time. D is the diffusion coefficient 
of the micelles. fj,s is the Stokes mobility (the electrophoretic mobility divided by the elementary charge), n is the viscosity of 
Isopar M. as is the micelle Stokes diameter determined from Eq. (|13[) . Q is the total charge of one sign in the cell determined 
from peak integration, no is the volumetric concentration of charges micelles, in nano-mol/liter. 



the cell is devoid of mobile charges since bulk generation 
is a relatively slow process. The ability to control the 
bulk charge density makes this system a unique physical 
laboratory. 

The spike-peak structure is the robust feature of the 
experimental and model transient currents. The devel- 
opment of the peak signals the onset of cell polarization. 
The reverse current is nonmonotonic because even after 
the voltage turn off, the electrodes carry a significant 
amount of charge. Initially, the latter keeps a significant 
portion of the mobile charges inside the cell close to the 
electrodes. As a result, the internal charges are released 
into the bulk not all at once, but gradually. At high ex- 
ternal voltages the peak shape saturates. All of these 
features can be reproduced in a simple one-dimensional 
drift-diffusion model. Such a model has been formulated 
in dimensionless units and solved numerically for a range 
of the two model parameters. The dimensionless peak 
time increases logarithmically with the charge density 
and lies in the interval 0.10 ± 0.05 for a wide range of 
parameters. Measured in real seconds, the peak time 
scales quadratically with the cell gap and inversely with 
the micelle mobility. The scaling provides a simple way 
to estimate the micelle diffusivity as D « 0.1 d 2 /t p . Tak- 
ing into account the skewness of the peak allows a more 
accurate determination of the micelle parameters, which 
have been summarized in Table [I] The Stokes diameter 
of the micelles is estimated to be 12.2 ± 0.5 nm. From 
the temperature dependence of the total charge density, 
the activation energy of 0.126 eV has been inferred. 
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